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Abstract 

The dynamics of moving vortex lattice is considered in the framework of the time dependent 
Ginzburg - Landau equation neglecting effects of pinning. At high flux velocities the pinning 
dominated dynamics is expected to cross over into the interactions dominated dynamics for very 
clean materials recently studied experimentally. The stationary lattice structure and orientation 



depend on the flux flow velocity. For relatively velocities V < V c = 7#o/7, where 7 is inverse 

diffusion constant in time dependent Ginzburg - Landau equation, vortex lattice has a different 
orientation than for V > V c . The two orientations can be described as motion "in channels" and 
motion of "lines of vortices perpendicular to the direction of motion. Although we start from 
the lowest Landau level approximation, corrections to conductivity and the vortex lattice energy 
dissipation from higher Landau levels are systematically calculated and compared to a recent 
experiment. 

PACS numbers: 74.60.-w, 74.40. +k, 74.25.Ha, 74.25.Dw 
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I. INTRODUCTION 



The static Abrikosov flux lattice has been experimentally observed since sixties by great 
variety of techniques and lateral correlations have been clearly observed recently up to tens 
of thousands of lattice spacings [1]. The remarkable advances in decoration, small-angle 
neutron scattering and muon spin rotation techniques allowed recently direct glimpse into 
the structure of the moving Abrikosov vortex systems [2-5]. It shows that at small flux 
flow velocities vortices move in channels as predicted in [7]. When the flux flow velocity 
increases beyond the one corresponding to the critical current, one observes a relatively well 
correlated hexagonal lattice. The channels and the plastic flow at relatively low velocities 
are explained by influence of pinning on the basis of theoretical arguments [8] and confirmed 
by numerous simulations [8-12]. At high velocity of the moving lattice (corresponding to 
high electric field), the influence of disorder is expected to diminish and a "moving Bragg 
glass" appears [8, 13]. Indeed Bragg peaks roughly at positions of the hexagonal lattice were 
observed [5] recently. 

Since the theoretical prediction of the moving Bragg glass exhibiting transverse peak 
effect [13], a lot of effort has been put into the simulation of the high driving force phase of 
the moving vortex system [10-12]. In particular it was found [10] that as the driving force 
increases (or disorder decreases) the vortex lattice suddenly changes orientation for a period 
of time and then returns to a "regular" drift mode. The main emphasis in these studies 
mentioned above is still the effects of pinning on the moving lattice. 

Experiments at low (below 100G) magnetic field and slow flux moving velocity (of order 
fim/sec) showed that the orientation of the moving vortex lattice is tied to the direction 
of motion, namely, when nearly hexagonal lattice is observed, one always observes the ori- 
entation depicted on Fig. la, never the "rotated" one of the Fig. lb [4]. Here the effect 
of pinning cannot be ignored and plays an important role in the orientation of the vortex 
lattice. However the most recent small-angle neutron scattering and muon spin rotation 
experiment can probe the moving lattice at much higher velocities of order cm /sec or even 
higher. The results about the orientation of the moving lattice obtained in [5] seem to be 
different from the case at low magnetic field and slow flux moving velocity. 

The effect of pinning is expected to be smaller at higher velocities. Alternatively one can 
ask what happens in very clean materials. A recent experiment in Pn — In seems to belong 
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to this category [5]. As the pinning influence diminishes with increasing flux velocity, it 
is natural to ask what happen in the limit of the highest possible flux velocity (of course, 
eventually the electric field destroys superconductivity, so that the mathematical limit of 
the infinite driving force is unphysical) disregarding pinning altogether. 

The question of the orientation of the vortex lattice usually does not arise in the static 
case. Without external electric field singling out a particular direction one has a complete 
degeneracy of possible orientations of the hexagonal vortex lattice. This is not surprising 
for a sufficiently symmetric material (like NbSe2 frequently used in experiments belongs to 
this category): the rotational symmetry ensures that the free energy is independent of the 
hexagonal lattice orientation. The rotational symmetry is broken by the motion of fluxons 
as was confirmed experimentally [4, 6]. Naturally one could ask whether the particular 
lattice orientation observed for example in [4] is necessarily tied to pinning or might appear 
in clean superconductors as well. Furthermore, the lattice also can be deformed though the 
deformation apparently is very small (see Fig.lc,d in [4]). Is there a deformation even before 
pinning centers disorder the lattice? 

It would be difficult to address the question of the moving vortex lattice structure using 
phenomenological models like the elastic medium [13] (in which individual vortices are simply 
not "seen") or approximating vortices in the London approximation by interacting lines or 
points Ti in 2D [12]. To give an example of the problems in the London limit, let us consider 
equations of motion for vortices. The driving force F is the Lorentz force and the dynamics 
is assumed overdamped: 



where U (r, — r 3 -) is the inter - vortex repulsive potential. The solution of these equations 
in the absence of pinning is obvious: the "boosted" hexagonal lattice of any orientation 
irrespective of the direction of F. Thus the orientation of the lattice depends solely on 
initial conditions, at least in the clean case. Therefore the approximations made in the 
above phenomenological approaches are too strong. 

In this paper we use the time dependent Ginzburg - Landau (TDGL) model to study 
the vortex motion and structure. The TDGL approach has been remarkably successful in 
describing various thermodynamical and transport properties [14]. Progress in obtaining the 
theoretical results from the model can be achieved only when certain additional assumptions 
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are made. One of the often made additional assumption is that only the lowest Landau level 
(LLL) significantly contributes to physical quantities of interest. The LLL approximation is 



moving lattice were performed at field far below the static H C 2(T), it has been shown long 
time ago [16, 17] that in the presence of electric field E the effective H C 2(T, E) = H C 2(T) — 
^ 2 \/ 2 $ /(87r) where V = cE/B is the velocity of fluxons and 7 is the inverse diffusion 
constant setting the time scale in TDGL approach. This field H c2 (T,E) could be much 
smaller at not very small fluxon velocities (electric field suppresses superconductivity even 
more effectively than the magnetic field). Therefore effectively one can move into the region 
of validity of the LLL approximation at sufficiently large currents. Moreover one expects 
that, even beyond the region of validity of the LLL approximation, physics is qualitatively 
the same. 

We solve TDGL equations for a moving vortex solid without disorder and find the vortex 
structure to which the moving lattice relaxes irrespective of initial conditions [16-18]. It 
turns out that the preferred lattice is rhombic. The distortion is velocity dependent. Re- 
markably the orientation is the same as on Fig. la, namely agrees with experiments only at 
velocities exceeding the critical one (of order of cm /sec for superconducting type II "low T c 
metals). Below it the orientation is rotated by 30°. 

The paper is organized as follows. Model is described, symmetries analyzed and perturba- 
tive mean field solution developed in section II. The general formalism is developed to treat 
the non Hermitian part of the equation. The shape and the orientation of the vortex lattice 
and the reorientation transition are described in section III. Then in section IV we calculate 
corrections due to higher Landau levels and derive general expression for conductivity. It is 
compared with a recent experiment. Section V is a summary. 

II. MODEL AND ITS PERTURB ATIVE FLUX FLOW SOLUTION 
A. Time dependent GL Model 

Our starting point is the TDGL equation [19] 



valid for H > 



H c2 (T) 



in the static limit [15]. Although most of the experiments concerning 
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The static GL free energy is: 



where a and b' are phenomenological parameters, 7 is the inverse diffusion constant which 
controls the scale of dynamical processes via dissipation. As usual the magnetic induction 



where the first term is the contribution of the normal liquid in the framework of the two 
liquid model and the second term is the supercurrent 



Tensor a n is the normal state conductivity. We assume that the coefficient of the covariant 
time derivative term 7 in eq.(2) is real although a small imaginary (Hall) part is always 
present [18]. The general case will be discussed in section V. 

We make several approximations (identical to those made in [20] and major parts of 
[17]) so that the problem becomes manageable. The physical conditions allowing those 
approximations are the following. Temperatures and magnetic fields are close "enough" to 
H c2 (T). Under this assumption the order parameter ip is suppressed compared to its Meissner 
value. In this paper we will also assume strongly type II superconductivity k — A/£ >> 1 
(£ 2 = h 2 / (2m a t,aT c ) , A 2 = ^JX^t )• Magnetic field is very homogeneous since the vortices 
overlap. Characteristic length describing the inhomogeneity of the electric field was identified 
in [17]: C 2 — ^p-A 2 and since typically a n ~ 7, thus ( » £ and the electric field is 
assumed homogeneous. Therefore the Maxwell type equations for electromagnetic field are 
not considered. The time independent vector potential will be taken in Landau gauge 
A = (By, 0, 0) and describes a nonfluctuating magnetic field in the direction — z. The scalar 
potential is also independent of time Aq = Ey and describes the electric field oriented along 
negative y axis. The vortices are therefore moving along the x direction. We neglect thermal 
fluctuations and disorder on the mesoscopic scale. 

Throughout most of the paper we will use the following physical units. Unit of length 
is the coherence length £, unit of magnetic field is H c2 = A = 4 " a ^'" , and the unit 
of energy (temperature) is T c . In these units the magnetic field is denoted by b = B / H c2 . 





(5) 
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The asymmetry of masses between the z direction and the x — y plane can be removed by 
rescaling coordinates and time: x — > £#/ \fh, y — > ^y/ Vb, z — > C,z/ ^Jbmjmab, t — > ^-t. The 
TDGL equations, after the order parameter field is rescaled as well ip — > ^y^ijj, is: 

= L^ + m 2 , (6) 

1 

2 



L = D t -l-[Dl + d 2 y +d 2 z ]-a, 



where a = 1 ~ T JJ' C , f = ^2Wy^^^ s scaled vortex velocity (in units of 2a/27i\B/$ /7) and 
covariant derivatives are defined by -D^ = J- — iy and D t — + ivy. Since 9^ commutes 
with L, the equations are invariant under the z translations, the z dependence of the solutions 
decouples and is generally a plane wave. It is easy to see that the relevant solution does 
not break this symmetry and is therefore constant with respect to z. Consequently we 
consider the problem as a 2+1 dimensional one (note however that if the 3D disorder or 
thermal fluctuations are included one can not ignore the z coordinate as the configuration 
of disorder can destroy the translational symmetry along the z direction) . 



B. Expansion of a nontrivial solution around dynamical phase transition point 

The line in parameter space (a,v), which separates the normal region in which the only 
solution is ip = from the flux flow nontrivial solution region, has been found by Hu and 
Thompson [17]. We will construct a perturbative solution of the TDGL equations near the 
mixed state - normal phase transition line analogous to the one in statics [21]. The range 
of applicability and precision of the LLL approximation at large k in statics was explored 
recently [15] . The main difficulty in the dynamical case is that the linear part of the equation 
L is not Hermitian due to the dissipation term D t . 

General idea of the expansion around a bifurcation point of a nonlinear equation is as 
follows. One looks for a set of eigenfunctions of the linear part of eq.( 6): 

The operator L consists of two parts: the usual Hermitian Hamiltonian of particle in mag- 
netic field — \ \D\ + d 2 ] and the anti - hermitian covariant time derivative D t . The complete 
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set of eigenf unctions with "quantum" numbers N and p x = p is: 

1 



exp[i(px — out)]H]sr(y — p + iv) exp 
2 



-\ {y-p + ivf 



1 v 

—a + N+- + — — i(u — vp) , 



(8) 



where Hn are Hermit polynomials. Unlike the usual case of a Hermitian operator, eigen- 
functions and eigenvalues of the Hermitian conjugate of the operator are different: 



Npuj 



e 



NpuiH>Npuj 
1 



exp [—i (px — ujt)]HN(y — p + iv) exp 



e 



Npu) 



1 V 2 

-a + N+ - + — + i(uj — vp) . 



(9) 



Note that <p is not a complex conjugate of <f). The orthogonality relations in the dynamical 



case involve both (f) Npw and 4> NpuJ : 



I ^ATpJ^y^jvw^'^'O = (2n) 2 5 NN/ 5(p-p')5(u-uj r ). 

J x,y,t 



(10) 



Npw (x,y,t)(f) Np0J (x,y,t)) xyt = 1, 



where the averaging over space and time is denoted by (■■■) xyt - 

The bifurcation (in this case the dynamical transition) occurs when there exists a set of 
eigenfunctions of L' with zero eigenvalues Qn P lu = 0: 



1 v 2 

a bif (v) = ^+2 + 7' 
us = vp. 



(11) 
(12) 



It is clear that solutions with N > are unstable as in the static case [19]. Equation (11) 
with N = gives the phase transition line of [17], while eq.(12) selects the "zero manifold" 
in the space of functions. We define the "distance" from the transition line 



a h (v) = a - a Mf (v) = a ~\~\ 



(13) 



When dh(v ) > 0, the nonlinear TDGL equation 



Lip + = L sh tp - a h (v)ifj + 



(14) 



L sh = L + a h (v) 



is solved perturbatively in with a nonanalytic prefactor, as in the static case: 



$ = [a h (v)} 1/2 [$° + o fc $ 1 + ...]. 



(15) 



To order [a^] , the equation linearizes: 



0. 



(16) 



Therefore $ belongs to the "zero manifold" and thereby can be expanded : 



N=0,p,w=vp — 




P' 



(17) 



p 



p 



with coefficients c p determined by the next order equation. As a result, since all the <f> p (x, y, t) 
depend only on the combination px — ut = p(x — vi) rather than separately on x and t, 
vortices move in the direction perpendicular to both electric and magnetic field with constant 
velocity v. To order [a?i] 3//2 , one obtains 



Multiplying this equation by <p p and integrating over (x, y, t) using of the orthogonality 
relation eq.(10) one obtains the following infinite set of nonlinear algebraic equations: 



We will study the solution of this set in the next section. 

III. SHAPE AND ORIENTATION OF THE MOVING LATTICE. 
A. Symmetry and energetics considerations 

It is well known in the static case that there is a solution of GL equations for any lattice 
symmetry. The same is true in the dynamical case as well, but the symmetries should 
take into account the motion of vortices. Define the covariant derivatives in a matrix 2+1 
dimensional form (summation over repeated indices assumed): 



L sh $ 




(18) 



Pi,P2,r 




(19) 
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and the Landau gauge 






1 

















— V 






(20) 



is used in our paper. All indices run over space /i — 1 (x), 2(y) and 3{t). The electromagnetic 
translation operators satisfying [7^, D^] = are: 

1 



JdP 



cxp 



3 idp 



(21) 



where generators are P M = —i (<9 M — ib u ^x u ) (note transpose in the matrix 6^). Operators 
p M = — i<9 M are usual (not "electromagnetic") translation operators. The following commu- 
tation relations 

[P M , P„] = i (V " M • (22) 

can be verified. Thus we will have [id r P, id 2 -P] = ~id la d 2 f3 {b a p — bp a ). Using the Hauss- 
dorf formula one checks that the electromagnetic translation operators obey e ld i p e * d 2 p _ 
e id 2 Pgidi P e [jdi P, id 2 P] jf { j i an j ^ are lattice vectors which preserve the symmetry 
of the system (when one translates system by di or d 2 , the system will be unchanged), one 



shall require e ldl ' p ifj = e^'^ip and it will lead to 



id 2 -P„ 



gidrPgUfc.P^, = ^.P^-Pg^.P, ld2 .P]^ = e [ idl .P, id 2 .P]^ = ^ 



Therefore we should demand 

[idi-P, id 2 -P] = i2ix x integer. 
This requirement is satisfied by the following basic translation symmetry vectors 

d«= ^(1,0,-1) 

d<°> = t(v,0,1). 



(23) 



Here oa is the lattice spacing along the direction of motion, r is arbitrary (a continuous 
translational symmetry). The flux quantization (one flux quantum per unit cell assumed) 
determines r'\ r'a\ = 2n. The d^ translation symmetry leads to discrete parameter 

2n l- l 

p= — i = gi 

OA 



in eq.(17), and the set of equations eq.(18) will take a form 



(24) 



exp |-i [(5-/1 + ivf + (gl 2 + ivf - v 2 ] j 
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It can be solved as in the static case by an Ansatz: 

ci = 

with the Abrikosov function 

9 



iirrl(l+l) 



Pa(v) 



exp {2%iriil2} 



x 



(25) 



11,12 



exp |-i [(5-/1 + ivf + (gl 2 + ivf - v 2 ] 



Consequently 



where 



$°(x,y,z) 



<p(x,y), 



-^(y - gl-iv)' 



(26) 



(27) 



v(x,y) = \ ^y^ ex P [ il (g( x - vt ) - nr ( l + i))]exp 
V V 

is normalized by (\<^\ 2 ) xy = 1- 

In the static case a solution which has minimal free energy is physically realized. The 
free energy is proportional to — [a h (0)] 2 /(2j3 A (0)) which therefore should be minimized. 
This means that one should minimize (3 A (0) . The minimal (3 A (0) = 1.16 is obtained for 
the hexagonal lattice. Similar reasoning cannot be applied to the moving lattice solution of 
the TDGL equation since the friction force is non conservative. Under these circumstances 
Ketterson and Song [22] calculated the work made by the friction force: 

27<|A^| 2 )^. (28) 



S = ^S 
dt 



The preferred lattice structure in the steady state corresponds to a state with largest S. For 
the lattice solution of TDGL equation one obtains to leading order in a h : 

S OC - ; ; X 



Pa(v) 



d 



v 2 \a h (v)\^ v 2 



I — + ivy j exp [il(g(x — vt) — nr(l + 1))] exp 



-\{y-9i-™f 



(29) 



2(5 A {v) 
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We therefore shall minimize (3 A as function of r and oa- This is consistent with the static 
case. 

B. The stationary orientation of the flux lattice. The reorientation transition at 
high flux flow velocity 

We found that the minimum of f3 A {v) appears always when r = 1/2, namely for rhombic 
lattices. Therefore from now on we consider these lattices only. As a function of an angle 
of the rhombic lattice tan# = 4n/a 2 A (see Fig. 1 for definition of 9) it generally has two 
minima, see Fig. 3. In the static case the two minima are degenerate with 9 = 60°, 
30° corresponding to perpendicular orientations of the hexagonal lattice, while for nonzero 
velocity the degeneracy is lifted. Note that originally [16, 17] it was assumed that the lattice 
is hexagonal also in the dynamical case. Generally the shape is not strongly distorted for 
physically realizable velocities. For velocities smaller than v c = 0.95 angle 9 close to 60° 
(the orientation of Fig. lb) is preferred over the one close to 30° (the orientation of Fig. la), 
see Fig. 3c. The dependence of the angle 9 on velocity can be very well fitted in the whole 
range v < 0.5 by 

9 = 30 - OAv - 2Av 2 . (30) 
The Abrikosov function also depends on velocity increasing according to 

P A (v)=f3 A (0)(l + l.25v 2 ), (31) 

where f3 A (0) = 1.1596 is the static value for hexagonal lattice. As the critical velocity is 
approached the two minima coincide, see Fig. 3b. Beyond that point the preferred structure 
is just the opposite, Fig.3a. The transition is first order and the coexistence region should 
exist. 

We now make A few comments about the orientation of the lattice. The reader might have 
noticed that the orientation of the lattice is not completely arbitrary since direction of the 
vector di coincides with the direction of the vortex motion. The most general f3 A {v) is given 
by eq.(25) with arbitrary r. We minimized numerically the Abrikosov (3 function and found 
that the solution with the largest dissipation is always of the more symmetric type r = 1/2. 
One can argue that despite the fact that electric field breaks the continuous rotational 
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symmetry, it still preserves a discrete transformation y — > —y,ip — > ip*. The solution 
r = 1/2 preserves this discrete symmetry. This symmetry is unlikely to be spontaneously 
broken. Indeed the symmetry was observed in the experiments (for example, in [4]). 



IV. NONLINEAR CONDUCTIVITY AND BREAKDOWN OF THE LLL SCAL- 
ING IN TRANSPORT 

In this section we first calculate the leading higher Landau level corrections to the solution 
of the TDGL equation eq.(6). Then we use it to derive the correction to the LLL scaling of 
conductivity [18, 20, 23]. 

A. Higher orders in correction to the moving lattice solution 

Using the same symmetry arguments as for the leading order, the second term in eq.(15) 
can be expanded as: 



N=0 



(32) 



lfN = y T^iV! ^2 ex P^ l (9(x-vt) -nr(l+ 1))] 
xH^(y — gl — iv) exp 
Multiplying eq.(14) by Tp N for iV > 0, one obtains: 

NC}, = -(3~ 3/2 (lp N <p*<p<p) (33) 
To find Cq we need in addition also the order cC^ 2 equation: 

L s/l $ 2 = $ 1 -(2$ 1 |$Y + $ 1 *$ $°) (34) 

Inner product with gives: 

oo 
N=l 

Note that for hexagonal lattice (Jp N (p*(p(p) ^ only when N = 6j, where j is an integer. 
This is due to hexagonal symmetry of the vortex lattice [21]. In statics (3 N = (Jp N ip*(pip) 
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decreases very fast with j: (3 6 = —.2787, (3 l2 = .0249 [15]. Because of this the coefficient of 
the next to leading order is very small (also an additional factor of 6 in the denominator 
helps the convergency) . 



B. The LLL scaling in nonlinear conductivity 



In the flux flow regime, in addition to the normal state conductivity, there is a large 
contribution from the Cooper pairs represented by the order parameter field. It was noted 
in [18, 20, 23] that the LLL contribution to nonlinear conductivity 

is proportional to the superfluid density. The scaled dimensionless conductivity is defined 
as a scaled = ^r-o" and a sca i ed in LLL approximation is: 

c *y 

a LLL = ^ (^iLLdy^LLL ' ^LLLdy^LL) = <|* ' LLL?) ■ (37) 

The last equality is due to the general property of the LLL functions, see eq.(8). It follows 
naive expectation of "Drude" like formula [19] with |\I/lll| 2 playing a role of "charge carriers" 
density (meaning here Cooper pairs). 

To leading order in using the results of the section II one gets 

2(3 A {v)v P A {v) 
where a,h(v) = (1 — t GL — b — v 2 )/{2b). At finite v there is an exponential factor coming 
from the nonorthogonality of eigenfunctions of a non Hermitian operator and, in addition, 
similar dependence in f3 A and quadratic in a^. In the limit v — > one recovers the Ohmic 
expression (see [18]) returning to standard units 

(l) l-t GL -b _ c 2 7 

aLLL = a ° WAV ' " = w' (39) 

while the leading nonlinear (cubic) correction is, using eq.(31), 

(3) tGL + b — 5 2 



*-=*° J feor'- (40) 



where v = £f 
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C. Leading correction to the LLL scaling 



Generally to all orders in ah one can write \1/ = J2 N Cn^^^n 

o= a ^^C* N (a h )C M (a h ) (tp* N d y (p M - ip M d y (p* N ) = vo^2C* N {a h )C M {ah)(TNM (41) 

V NM NM 

For N > M and N — M even integer 



»™ = " \/ 2 ^(-" 2 )^[- 2 iK ,+1 (-2» 2 ) + ^iK'-'M^K", (42) 



where L(y) are Laguerre polynomials. This contribution is always subohmic (?nm ~ v N ~ M 
at small i>. If TV — M is odd, the contribution vanishes. The diagonal contributions are 
simpler 

a NN = [Lj v _ 1 (-2t; 2 ) + Lj v (-2^; 2 )]e ^,2 (43) 

and have an ohmic part 

a NN = 2N + 1. 

The first term, proportional to Landau orbital number N, is responsible for the breaking 
of the naive "Drude" like expectation that conductivity is proportional to |\I/| 2 [19]. One 
observes that higher Landau levels contribute to conductivity more than to . One can 
interpret this as "increased charge movers density". 
Thus the ohmic conductivity has two contributions: 



C7« 



a ^(2iV + 1) \C N (a h )\ 2 = a 1 + a 2 ; (44) 

N 

a 1 = a (\q\ 2 ); a 2 = 2 a ^ N \C N {a h ) | 2 , (45) 

N=l 

the first proportional to the superfluid density, while the second, the HLL part, is not and 
is of order a\ only. Substituting expressions for C from the previous section, we obtain for 
the Ohmic conductivity to order a\ 



<J]_ = <7 



a h „ o|_ \^ 0N_ 



(46) 



where all the quantities are taken in the limit v — > 0. The sum rapidly converges in the 
static (or low velocity) case: J2n n = 0-0131. 
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D. Comparison with experiment 



On Fig. 4 we compare the results with a recent experiments at high currents (electric 
fields) of [24] on Nb in which vortex velocities as high as 10 5 cm/ sec. We used the same 
values of the Ginzburg - Landau parameter k — 9.4 and the inverse diffusion constant 
7 = 1.17 sec/cm 2 to fit all three curves corresponding to magnetic fields H = 80 mT, 100 
mT and 120 mT for "cold" sample with T c = 8.6 K. We used the measured (inset on Fig. 2 
of [24]) H c2 = T c dHc ^ T) \ T=Tc = 4.4 T. The temperature was T = 7.8 K close enough to 
T c so that the a\ correction was always below 10%. The value of parameter 7 is in good 
agreement the measured normal state resistivity of 9.9 fiQ cm. The results agree well with 
the flux flow Ohmic conductivity data at relatively low currents (still well above the critical 
current) exhibiting the 1/H behavior presented in Fig. 2 of [24]. 

One observes that the full expression (solid lines) is closer to the experiment at very high 
electric fields. Several curves for magnetic field ... are given. The smallest is clearly off the 
LLL approach range. 



V. CONCLUSION 



To summarize, we have considered the dynamics of vortex lattice neglecting effects of 
pinning. We studied the time dependent Ginzburg - Landau equation in the lowest Lan- 
dau level approximation. The validity region of the LLL approximation, as in the static 
case, we require ah = jb H c2 (T) — B — c << 6, the factor 6 coming from cancel- 

lations of the higher Landau level effects due to hexagonal symmetry (even the hexago- 
nal symmetry is approximate in the moving lattice). We systematically calculated higher 
Landau level corrections to conductivity and the vortex lattice energy dissipation. The 
stationary lattice structure depends on the flux flow velocity. While for small velocities 
V < 2f C A/27ri?/$o/7 ? v c = 0.95 vortex lattice is oriented as in Fig. lb, while beyond this 
velocity orients like in Fig. la. We emphasize that in our calculation the pinning effect was 
disregarded. Of course, as was firmly established in numerous theoretical and experimental 
investigations, pinning significantly can modify the picture for low velocities. Pinning gen- 
erally "prefers" configuration of Fig. la and this is a possible reason why the experimental 
observed orientation is depicted as Fig. la. However one can expect that for higher velocities 
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and very clean samples the pinning dominated dynamics crosses over into the interactions 
dominated dynamics considered in the present work. The high velocity of order cm/ sec is 
unlikely to be seen in decoration experiments. However other techniques like SANS and 
muon spin rotation [5] and possibly Lorentz microscopy [25] are able to detect the lattice 
structure even at such relatively high velocities. At very high velocities the results for 
nonlinear conductivity agree with recent experiments [24]. 
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Figure captions 

Fig. 1 

Two possible orientations of the (approximately) hexagonal vortex lattice, (a) direction 
of the flux lines is the same as the nearest neighbors lattice orientation, (b) direction of the 
flux lines is perpendicular to the nearest neighbors lattice orientation. 

Fig.2 

The flux lattice geometry: , d( 2 ) are the translational symmetry vectors which deter- 
mines the primitive cell of the flux lattice. The angle between these two vectors is 9. 

Fig.3 

Dependence of the Abrikosov f3 parameter on orientation and shape of the vortex lattice 
moving with scaled velocities v = 0.5.0.95, 1.1. The angle 9 is defined as an angle between the 
direction of motion and a crystallographic axis in direction of the symmetry transformation 
d 2 . The minimum favors the smaller angle close to 30° corresponding to structure of the 
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Fig. la for v < v c , while the other local minimum corresponding to Figlb (angles close to 
60°) is preferred for v > v c . 

Fig.4 

Current - voltage curves at high flux flow velocities. Data of ref. [24] on Nb films at 
T = 7.8 K (symbols represent different magnetic fields) are compared with theory combining 
the linear (Ohmic) contribution eq.(46) and the cubic correction eq.(40). 
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Li, Malkin and Rosenstein, Fig.1 a 




Li, Malkin and Rosenstein, Fig.1 b 




Li, Malkin and Rosenstein, Fig. 2 




Li, Malkin and Rosenstein, Fig.3 




Li, Malkin and Rosenstein, Fig. 4 




